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Abstract 

In the context of a dynamical Ginzburg-Landau model it is shown numerically that 
under the influence of a homogeneous external current J the vortex drifts against the 
current with velocity V = —J in agreement to earlier analytical predictions. In the 
presence of dissipation the vortex undergoes skew deflection at an angle 90° < 5 < 180° 
with respect to the external current. It is shown analytically and verified numerically 
that the angle 5 and the speed of the vortex are linked through a simple mathematical 
relation. 
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1 Introduction 



In a recent publication it was proposed to study the dynamics of isolated Abrikosov 
vortices in the framework of a phenomenological effective field theoretic model , which 
alternatively can be viewed as a time dependent version of the Ginzburg-Landau equation 
(TDGL). Therein, it was possible to derive analytically the equation of motion of the 
guiding center of the vortex (Hall equation), under the influence of any kind of external 
force 0. According to the latter, the guiding center of the vortex moves in a direction 
perpendicular to the externally applied force in analogy to the planar motion of charged 
particles in the presence of a perpendicular magnetic field. The importance of the Hall 
equation stems from theoretical considerations which suggest that the guiding center 
describes quite accurately the "mean" position of the vortex, a suggestion which was 
verified and quantified by a recent numerical study Q]. It is thus possible to obtain 
useful information concerning the motion of a vortex by simply invoking Hall equation. 
For instance, by plain implementation of the latter one finds that in the presence of a 
homogeneous external current the vortex drifts against the current. 

By construction, the Hall equation is limited to the description of the gross features 
of the motion of the vortex. On the other hand, an understanding of the finer details 
of its motion requires a detailed solution of the TDGL equations. Such an undertaking 
has to rely on numerical methods, due to the nonlinear nature of the relevant equations. 
The main purpose of this paper is to report on a numerical study of the dynamics of 
vortices within this model. In particular we will study the response of a vortex to an 
external current with or without dissipation, an issue of obvious interest in the context of 
superconductivity 

The paper is organized as follows. Section 2 contains a general introduction to the 
model. Its relevance to the physics of the superconductor is commented upon, and the 
main theoretical predictions concerning the motion of the vortex are illustrated. Finally, 
it is shown that an alternative interpretation of Magnus force 0] arises naturally in terms 
of the Hall equation. In Section 3 we incorporate the effect of dissipation and that of an 
external electric current in our field theoretical formalism and we derive the equation of 
motion for the vortex through explicit calculations. The results of an extensive numerical 
study are presented in Section 4. Vortices are found to behave in accordance to our 
earlier theoretical predictions. As a byproduct of this study we investigate the detail 
of the internal "Cyclotron motion" of vortices and we show that an isolated vortex is 
spontaneously pinned. In the concluding Section 5 we discuss the relevance of our results 
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to the actual superconductor and we propose more realistic 3-dimensional studies. 

2 The Model 

Most of our techniques and conventions are described elsewhere so we briefly outline 

here the physical and mathematical background of the theory. The model in question was 
originally introduced by Feynman as a natural dynamical extension of the Ginzburg- 
Landau static theory of superconductivity. Indeed, by attributing the correct physical 
content to flelds and parameters, it becomes a rather realistic phenomenological model of 
a superconductor ||^. The model admits inflnitely long, smooth, cylindrically symmetric 
flux vortex solutions, whose static properties together with the properties of pairs of them, 
have been studied in detail in reference Our objective is to study some aspects of 
vortex dynamics within this model. In doing so, we will ignore excitations along the axis 
of the vortex (which by convention is taken to be parallel to z-axis) and we will directly 
deflne the model in two space dimensions. 

As usual, the important dynamical variable is a complex order parameter \& which 
may be thought as an electrically charged fleld coupled to the electromagnetic potential 
{Aq, A). The flelds satisfy the coupled system of TDGL equations (to simplify notation, 
flelds and coordinates are rescaled along the lines of reference @]) 

i4f = -^D^^ + + ^K^(^*^ - 1)^ 

^Ei = ei,d,B-J^ = p (2.1) 

with B = eijdiAj, Ei = —dtAi — diAo, Di = di — iAi, P, k are free parameters, and the 
spatial indices i,j range from 1 to 2. The supercurrent density J'* and the charge density 
p are given by J/ = ^[\E'*Dj\E' — c.c] and p = \E^*^ — 1. Note that in order to allow the 
possibility of a condensate (|\I'| = 1) at inflnity we introduce a background (positive- ion) 
charge density pb to neutralize the system. For simplicity is taken to be constant and 
homogeneous throughout i.e. p^ = —ps where ps is deflned as ps = \E'*\I'|oo = 1- 

The nonlinear system of equations (|2.1|) admits static, axially symmetric localized 
vortex solutions. They are similar in nature to the well known Abrikosov vortices of the 
static GL theory. Moreover, in a certain region of the parameter space their characteristic 
lengths (penetration depth, coherence length) fall well within the scale of a typical type II 



3 



superconductor However, the vortices in this model differ from the ordinary Abrikosov 
vortices mainly in one respect; although they carry zero electric charge as a whole, there 
is a local charge modulation in their interior i.e. they have non-vanishing charge density 
and as a result an electric field is associated with them. 

A remarkable feature of the model in hand, is the unusual response of the vortices to 
external probes. Indeed, as it was shown in references the motion of the vortex is 

governed by the equation 

Vii = ^ii Fi (2.2) 

where is the mean velocity of the vortex and F is the sum of the external forces 
acting upon it. N is the well-known integer-valued winding number or topological charge 
characterizing any finite energy configuration which counts the number of times the phase 
of \1/ rotates around the internal circle as we scan the circle at spatial infinity P] . is a 
conserved quantity and can be written as the integral of a properly chosen "topological 
density" r. Among other possibilities r may be defined as r = ^B, a definition which 
entails the familiar magnetic flux quantization and it is of obvious physical interest while, 
for our purposes the most useful form of r is 

r = [ekiiDk^YiDi^) - tB{^*^ - 1)] (2.3) 

The importance of the above formula stems from the fact that it appears in the expressions 
for the momentum and the angular momentum of the theory. Indeed, it was pointed out 
in reference that the Noether expression for the linear momentum of the model is 
ambiguous for any configuration with non-zero topological charge. The unambiguous 
expression turns out to be 

a, = .../A(2...r + i£.B) (2.4) 

which in turn leads to a radical revision of the physical interpretation of the momentum. 
The presence of the first moment of the topological density in the expression for the linear 
momentum, directly associates the latter with the position of the vortex. In fact, it is 
possible to show that the "mean position" of the vortex is described by a quantity 
R called the guiding center of the configuration and defined as Ri = —^2^^ijPj ■ By 
identifying R with the position and R with the mean velocity of the vortex, and 
assuming that a generic force Fj = ^ acts on the system, we end up with equation 
12T 



Contrary to intuition and Newtonian Mechanics, equation (|2.2|) implies that a vortex 
moves at a constant calculable speed in a direction perpendicular to the applied force. 
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This kind of behavior is analogous to the planar motion of charge particles under the 
action of a perpendicular magnetic field |^ and thus we call it Hall behavior and from 
now on we will refer to eq. (|2.2| ) as Hall equation. Vortex dynamical behavior becomes less 
exotic and Newtonian Mechanics reestablishes by adopting an alternate interpretation of 
the Hall equation. To do so, we rewrite eq. (|2.2| ) in the form : = — 2~]v ^ where 
Vi, is the velocity of the vortex in the x — y plane and F is the total force acting upon 
it. Multiplying both sides of the former equation by 27riV and taking the cross product 
by we get 

27rA^VL X = F (2.5) 



According to eq.( |2.5| ), for a vortex which moves with constant velocity Newton's Law 
X] F = is restored if we assume that in addition to any other force acting on the vortex 
another 'new' transverse force F^^ = —2'kN\'i x also acts on it. Bearing in mind that 
the total magnetic flux 0o equals 2ttN we see that Fj^ has similar form to the familiar so 
called 'Magnus force' |0| which is usually invoked to describe the motion of vortices in the 
superconductor. The similarity becomes more obvious in full units where this additional 
force reads Ft = — Ps^oV/, x e^- 

As we have explained in some previous work Hall behavior (and consequently the 
transverse force Fy ) is a generic characteristic of soliton dynamics in systems with non 
trivial topology and spontaneously broken Galilean invariance, due for instance to the 
presence of a crystal lattice, and has a clear mathematical origin. Yet one would like to 
have a more physical explanation for the appearance of F^ in eq. ( p. 51) . In fact, one can 
attain such an explanation by attributing the origin of Fy to the interaction between 
the magnetic flux of the vortex and the internal electric currents which are generated 



by the motion of the vortex [|10[. Specifically, let us assume that the vortex is moving 
with constant velocity V^, . In our field theory prescription a moving vortex with velocity 
\l is a field configuration \I'(r — VLt)e*^'^^''^^'^^'^*\ A(r — V^^t) (see ref. for more 
details on Galilean boosts in 2-d vortices). Let us now switch to the reference frame 
where the vortex is still. In that frame, the background ions of the crystal lattice form 
a homogeneous current of negative charge carriers with charge density —ps and velocity 
— Vi. This current interacts with the magnetic field of the vortex and as a result feels a 
Lorentz force acting upon it. Consequently the vortex feels a backreaction force opposite 
to the Lorentz force which can be easily computed and is found to be exactly the transverse 
force Ft mentioned above. 
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3 External Current and Dissipation 



To study the response of the vortex to an externally prescribed current J'^^*(x, t) we 
simply substitute J'' — > J* + J*^^* in eq. ( |2.1|) . Because of the external current the linear 
momentum (|2.4| ) of the system is no longer conserved. A straightforward application of 
the equations of motion yields 

j^Pk = Ft°''^'' = - j d^x ekiJr\x, t)B{x, t) (3.1) 

Assuming that the external current is uniform throughout the plane i.e. J^^*(5', t) = Jq, 
the only space dependent quantity on the right hand side of this equation is the magnetic 
field, and its integral, the total magnetic flux, is equal to 27cN. Thus, the equation of 
motion for the momentum takes the simpler form Pk = —27TNekiJoi- Correspondingly 
the time evolution of the guiding center or the "mean position" of the vortex reads 

which quite surprisingly implies that the vortex drifts against the current with constant 
calculable speed. Note that in a similar model, Manton |rT| arrives to the same result by 
considering Galilean boosts on vortex configurations. 

The pull of the vortex by the current has a simple explanation in the context of Hall 
equation. The magnetic flux of the vortex exerts on the electric current a Lorentz force 
in the —90° direction with respect to the current. Consequently the vortex feels a back- 
reaction force in the +90° direction and naively one would expect the vortex to move in a 
direction perpendicular to the current. However, as it follows from the Hall equation the 
vortex moves at +90° with respect to the applied force, and therefore, the vortex moves 
in a direction opposite to the external current. 

The introduction of dissipation in the system is a complicated task. Formally the effect 
of dissipation in the system is studied by adding a, phenomenological friction term in the 
TDGL equations. However, no such term has been derived on the basis of solid physical 
reasoning. Yet, one must have in mind that there are several restrictions in the form 
and the properties of any such term. Any friction term inserted in ( p.lD should meet 
the following conditions : a) it should vanish for any static vortex solution, b) it should 
decrease the total energy W of the system i.e. —dW/dt should be positive definite and 
finally c) should preserve the electromagnetic U{1) gauge invariance of the system, or 
equivalently it should preserve the continuity equation. Our choice, though not unique 
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was the most natural among a small set of candidates, and that was to add a term of the 
form Cd^ijdjB -with Cd a positive constant- on the right hand side of the equation of 
motion for the electric field in (12.11). 



By construction the friction term is gauge invariant and vanishes for any static configu- 
ration. The time derivative of the energy reads dW/dt = ^ / d'^x{Ei (Cd^ijdjB)} which 
by integration by parts becomes dW/dt = — ^ / d'^x{{eijdjEi) B} = — ^ / d'^xB'^. Thus 
we conclude that the friction term we propose meets all the restrictions mentioned above 
and therefore it seems to be a reasonable candidate for a phenomenological study of 
dissipation. 



:Pk = Cdj d^x {ekieimdmB{x,t)B{x,t)} = Cd j d^x {B{x,t)dkB{x,t)} (3.3) 



-Pk = -Su.Vlj with Sk, = Cdj d^x{dkBd,B} (3.4) 



The time evolution of the linear momentum is 
d_ 
dt' 

To proceed further we adopt the rather plausible assumption that in the presence of a 
suitably chosen external current as well as dissipation a steady state is eventually reached 
in which the vortex moves rigidly with constant velocity Vl. We assume further that the 
profile of the magnetic field associated with the steady state of the vortex is of the form 
B{x,t) = B{x -VLt^Vi). Then B{x,t) = -VLjdjB and eq. (pD reduces to 

d_ 

dt' 

A simpler version of this relation is obtained by invoking some further assumptions about 
the steady state profile of the vortex, which may be viewed as reasonable approxima- 
tions in the limit of a weak external current. Thus, we assume that the vortex retains 
approximately its initial shape and correspondingly its axial symmetry i.e. in the rest 
frame of the vortex the magnetic field is of the form B = B{p) , where p is the radial 
coordinate in the usual polar variables. Then the diagonal terms of the Skj tensor vanish 
while Su = S22 = Cd/2 J d^x {dpB'^) = r] and equation (|3l^) reduces to 

^Pfc = ~r]VLk (3.5) 

with 1] being a positive constant number. Eq. ( p.5|) implies that the effect of such a term 
on an axially symmetric vortex configuration moving with velocity is a force linear to 
the velocity of the vortex which opposes the motion of the latter. We thus believe that 
with the insertion of this term we correctly incorporate the effect of friction in our model. 
Finally, the equation of motion of the guiding center of the vortex in the presence of both 
a homogeneous external current and dissipation reads 

^Lk = = —Jok + ^^j^^^kmV^Lm (3.6) 
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Equation ( p.6|) implies that a vortex in the presence of a transport current, simply drifts 
against it (Vl = —Jo) ; this leads to a perfect Hall effect but with opposite sign to 
that of the normal-state [0. In the presence of dissipation the velocity of the vortex 
acquires a component in the perpendicular to the applied current direction. However, 
its longitudinal part still has opposite direction to the transport current, which in turn 
results in a sign change of the Hall effect [Q. By plain implementation of eq. ( |3.6| ) we 



find that the defiection angle 6 between Vl and Jo decreases from 180° to 90° as the 
contribution of the drag force increases. Furthermore, a simple relation for the deflection 
angle 6 is obtained, namely 

I I 

cos6 = — — - — - (3.7) 
I Jo I 

The former relation becomes obvious by simple inspection of flg.l. There are two forces 
acting upon the vortex ; i) the Lorentz force (ir^o^^nte _ 27TNekiJoi) and ii) the dissi- 
pation force (Frf = — r/V^). According to eq. (|2.2| ) the vector sum of and 'P Lorentz 
rotated by +90° and divided by 2t:N must be equal to the velocity of the vortex V^. 



The latter condition leads to relation (|3.7| ). 

From our analysis we conclude that within our model, and with or without dissipation 
the longitudinal part of the velocity of the vortex always has a direction opposite to the 
transport current which in turn results that the Hall effect in the vortex state will have 
a sign which is opposite to that of the normal state. We should note here that equation 
( |3.6| ) is not original at all. The major contribution of this work is its fleld theoretical 
derivation and its interpretation. Indeed, by multiplying both sides of (|3.6|) by 2TiN and 



taking the cross product with Cz we can write the latter in the form 

27riVVL X = -27rA^ Jq x - rjV l (3.8) 

All the terms above appear in the most common phenomenological theories of vortex 
motion [0| . The left hand side of (|3.8|) is the familiar though controversial Magnus force 



while the flrst term of the right hand side is the so-called Lorentz force [|1^ , or according 
to some other authors []T3| an inseparable part of the Magnus force. Finally, the last term 
of the right hand side of ( p.8|) is like the viscous drag force of the Bardeen-Stephen (BS) 



model |]T3[. However, one should notice that the Lorentz force comes with opposite sign in 
eq. (3^) with respect to what is commonly accepted in the literature, a difference which 
is of critical importance. In fact, the sign inversion of the longitudinal part of the vortex 
velocity in our analysis, has its origin in the sign change of the Lorentz force. The reason 
for the latter change is in the way we incorporate the external current in our model. We 
think of the current as a totally substantive object which interacts with the vortex. It is 
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like a vortex moving in a plane under the influence of a current which flows in a parallel 
plane just above the vortex plane, instead of a current flowing in the plane of the vortex 



and formed by the same carriers as those of the vortex |T3|. At first sight the second 
approach looks more natural. However, if we think of the vortex as a 3-dim object i.e. 
as a flux tube formed in a 3-dim superconducting film and interacting with supercurrent 
which is formed only on the top and the bottom of the film, our approach becomes more 
plausible and realistic. 



4 Numerical Results 

Our next objective is the numerical investigation of the dynamical behavior of the vor- 
tices in the presence of an external electric current and dissipation. Our computational 
techniques are described in detail elsewhere [Q, so here we present only a brief overview. 

To simulate the motion of a vortex, we first determined numerically the static profiles 
of the condensate and gauge fields characterizing a vortex of winding number 0. 
With these in hand, we laid down on the lattice a configuration of a vortex centered at 
(x, y) = (2, 0) at time t = 0. In order to maintain (as much as possible) the symmetries of 
the continuous system in its discretized form, we have resorted to techniques from lattice 
gauge theory [1^. The degrees of freedom were discretized on a spatial lattice so as 
to maintain exact (lattice) gauge invariance. The time evolution was implemented by a 
finite difference leapfrog method using equations ( p.l|) with the second one modified to 

^E, = e,,d,B - - + Cde,,d,B (4.1) 

in order to incorporate the effect of an external current and dissipation. The gauge 
freedom of the equations of motion was eliminated by imposing the temporal gauge Aq = 
0. The external current J*^^* was taken along the x-direction and uniform throughout 
the whole plane i.e. J^^* = Joe^. 

In any numerical calculation where partial derivatives are involved, the imposition 
of appropriate boundary conditions is a very delicate task. Here, the presence of an 
incoming and an outgoing external current at ± x-infinity, made this issue even more 
complex and forced us to use different boundary conditions at x and y-boundaries. At 
the y— boundaries of the film we imposed free boundary conditions by setting the covariant 
derivative in the normal to the boundary direction equal to zero {Dy'^ = 0). There are 
more than one ways to impose such a condition, and our choice was to set dy'^ = 
and Ay = 0. Also, in order to get vanishing magnetic field at the boundaries we set 
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dyArc = 0. At the x— boundaries we successfully imposed two different sets of b.c. There, 
in deriving the boundary conditions we took special care in order to preserve the discrete 
gauge invariance of the system, or equivalently the discrete version of the continuity 
equation, at the x— boundaries. Specifically, in analogy to the y— boundaries we set 
dxAy = 0, = and then, by imposing the continuity constraint we got 4[\E'*9a;\E'] = 
— Jo- As a consequence of the above b.c. the value of the x-component of the supercurrent 
at ±a;- infinity was fixed and equal to — Jo- Even though, this is what we expect to happen 
away from the vortex, still the b.c. for \E' sounds quite artificial and too constrained. 
Thus, we tried and finally used another more "natural" set of b.c. Namely we set d^A^ = 
dxAy = and the constraint of preserving the continuity equation yielded D^^^ = as 
a b.c. for \E'. Using both sets we got essentially identical results, and these have been 
exhaustively checked to be free of any boundary contributions. 

The simulations were done on a square lattice of 201 x 201 sites with a lattice spacing 
a = 0.15. The width (diameter) of the vortex was typically of the order of 2 in rescaled 
length units. The finite time step was chosen to be much smaller than the lattice spacing 
typically of the order 10~^ . To test our results we ran simulations in bigger lattices 401 x 
401 with the same or smaller lattice spacing, say a = 0.1, and the results obtained were 
all perfectly consistent. All our simulations were performed on various HP workstations at 
the University of Crete. A typical run of duration T ^ 100 time units, with At = 0.002 
on a 201 x 201 lattice needed about 15 hours of CPU time on a HP- 735 machine. 

Apart from the existence of the external current and the new boundary conditions, the 
algorithm we used here was identical to the one previously used in the numerical study of 
a vortex pair dynamics 0|. There, it turned out that the algorithm was extremely accu- 
rate. Here, as a sort of calibration of the algorithm, and mainly to avoid any systematic 
contributions from the b.c. we initially tested it in a system where the result is known 
analytically, namely at the trivial sector = 0. There, it is easy to see that the field con- 
figuration ^'(x, t) = exp[— I Jpt], Axi^, t) = Jq, Ay{-x., t) = = Ao(x, t) is a solution of the 
equations of motion in the presence of the external current with E(x, t) = = -B(x, t) and 
J**(x, t) = — J^^'*. We ran a preliminary simulation using as initial configuration the trivial 
vacuum = l,Ai = = Aq and turning on the external current at t = 0. Our naive ex- 
pectation was to see the system relaxing to the vacuum- current solution described above or 
to some gauge transform of it. Yet, after a short transient period the system dynamically 
relaxed to a time-dependent configuration where both the electric field and the supercur- 
rent oscillated vividly around their mean values < E(x) >= 0, < J^(x) >= — J^^'*. These 
oscillating modes could be attributed to the abrupt turning on of the external current. 
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To eliminate them, we introduced dissipation in the system and in a subsequent run we 
saw the system relaxing to the vacuum-current solution within a 10~^ accuracy. (Note 
that we introduced here a dissipation term of the form Ei = —C'^Ei + .... which does 
not meet all the criteria mentioned in the previous section and thus it is not appropriate 
in the vortex sector). The remarkably accurate convergence of the initial configuration 
to the vacuum-current solution provides a strong confirmation for the accuracy and the 
reliability of our algorithm. A byproduct of these runs is the conclusion that the response 
of the ground state to the application of an external electric current is the formation of a 
supercurrent which in average is equal and opposite to the latter. In fact, it is reasonable 
to assume that even in the N ^ sector, away from the vortex, the system will response 
in the same way. 

We switch now to the study of the vortex sector and we particularly consider the = 1 
sector. To study the dynamics of the vortices we carried out numerous simulations and 
the results confirmed with quite impressive accuracy the predictions of the theoretical 
analysis. We also experimented with the values of the parameters k, j3 and the dynamics 
of the vortices showed little sensitivity to those values. Indicatively we quote here the 
results of some simulations for the specific choice k = 2 and j3 = 1, which belong to a 
parameter regime that we believe to be appropriate for the description of the physics of 
type II superconductors 0]. Similar results though, were obtained for a large variety of 
the values of the parameters. 

We performed a series of numerical experiments for various values of the dissipation 
constant Cd and for a fixed value of the external current Jq = 0.025. The total duration 
of each simulation was 240 time units. The corresponding trajectories of the guiding 
center of the vortex are displayed in fig.2. As arises from the plot, in the absence of 
dissipation {Cd = 0), the guiding center of the vortex performs a rectilinear trajectory 
along the negative x-direction. Moreover, it is displaced by 6 space units (i.e. three times 
its diameter) in 240 time units, which corresponds to a mean velocity V/, = —JqCx- This 
behavior verifies the theoretical prediction ( p.2| ) qualitatively as well as quantitatively. 

To arrive to eq. ( p.2|) we assumed that the external current is homogeneous throughout 
the whole space. This assumption may take a weaker form; instead of a current occupying 
the whole space, we may introduce one which has a strip shape, i.e. it takes a non-zero 
constant value inside a strip of certain width, while it vanishes in the outside region. In 
principle, when the width of the strip is much larger than the size of the vortex, the vortex 
essentially realizes a homogeneous current all over space and responds accordingly. This 
assumption was tested and verified in our study. Specifically, we repeated our runs for 
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Cd = using an external current of the form 



r/' = o JT' = Jof{y) (4.2) 

where the function / is approximately equal to unity over a strip of width 2L and drops 
to zero very quickly outside that strip. A function which meets the above description is 

f{y) = e-d^l/^)" (4.3) 

especially for large n. The numerical calculations showed that any strip current with 
L > 6 has the same effect on the motion of the vortex to that of a homogeneous current 
occupying the whole space. Quite surprisingly, even for narrower strips, L < Q the 
trajectory of the vortex remains identical while its speed reduces. Indicatively, we found 
the ratio — ^ = 0.89, 0.95 and 0.98 for L = 3, 4 and 5 respectively. 

For non-zero values of Cd, after a small transient period the vortex relaxes to an almost 
rectilinear motion at an angle to the external current different from 180°. Indeed, the 
deflection angle 6 between the trajectory of the vortex and the direction of the current 
takes values from 180° to 125° with 6 decreasing as Cd increases. As arises from the 
lengths of the trajectories in fig.2, the measure of the velocity of the vortex is also a 
decreasing function of Cd in accordance to the physical interpretation of dissipation. 

An important issue to check is to which extent the motion of the vortex meets the 



assumptions we adopted in section III while deriving eq. ( |3.7| ). To see whether a steady 
state is eventually reached, we plot in fig. 3 the time evolution of the x and y-components 
of the drift velocity of the vortex for Cd = 4. There, we see that initially {t = 0), 
Vx = —0.025 = —Jo and V^^ = as if there was no dissipation at all. As the vortex 
moves on, dissipation turns on and its effect results in a non-zero transverse component 
of the vortex velocity. After a small transient period {t = — 4), the velocity of the 
vortex sets in an oscillating mode around a constant mean value. Strictly speaking the 
vortex does not seem to develop a steady state, but it is reasonable to assume that the 
contribution of the oscillating part is auto-canceled in average, and thus in a wider sense 
we can assert that the vortex finally relaxes to a steady state. Also, a detailed examination 
of successive level contours of the energy density establishes that the vortex moves quite 
coherently and retains its initial shape during the evolution of the simulations. It is thus 
quite interesting to question whether the trajectories shown in fig.2 satisfy eq. ( p.7| ). To 
calculate the defiection angle and the (mean) velocity Vl of the vortex we process the 
numerical data so as to linearize the slightly wavy trajectories of fig.2 and we cut out the 
initial part of the data which corresponds to the "transient period". After the relevant 
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calculations we tabulate the results in table I, where it is demonstrated that relation ( p.7|) 
is satisfied with quite impressive accuracy. 

Table I 





—cos5 


Vl 
Jo 


cosd 
Vl/Jo 


0. 


1.000 


1.000 


1.000 


1. 


0.984 


0.981 


1.003 


2. 


0.942 


0.934 


1.009 


4. 


0.816 


0.803 


1.016 


8. 


0.579 


0.582 


0.995 



It is of crucial importance for the present study to determine to which extent the de- 
tails of the motion of the vortex follow the motion of its guiding center. To investigate 
the latter we must somehow define the position of an extended object such as the vortex. 
Thinking of the vortex as an energy lump and given that this energy lump moves coher- 
ently retaining its initial axially symmetric shape, we believe that the maximum of the 
energy density (MED) of the vortex configuration accurately describes the real position 
of the vortex. Therefore, in our calculations we also keep a record of the location of 
the MED of the vortex configuration as time evolves. Note that there are at least two 
alternative definitions of the position or the center of the vortex, namely the maximum of 
the topological density or the position where \E' vanishes. All three of them yield similar 
results so here we present only data for the MED. The trajectories of the MED of the 
vortex (solid line) and of the guiding center (dashed line) of the vortex for Cd = and 
for various values of the intensity Jq of the external current are plotted in fig.4. We see 
that while the guiding center simply drifts in the negative x-direction, the motion of the 
MED is more complicated. In average it follows the motion of the guiding center, but 
its trajectory is modulated by an oscillatory pattern. This modulation is not a numerical 
effect but is an inherent characteristic of vortex dynamics. We have already encounter 
similar oscillating patterns in an earlier work while studying the dynamical evolution of a 
pair of vortices . These patterns are reminiscent of the motion that electrically charged 
particles perform in the presence of a magnetic field, which happens to be the prototype 
physical system which exhibits Hall behavior. Borrowing the terminology from the latter, 
we will refer to this finer motion with the name "cyclotron" . As in the original case, the 
amplitude of the cyclotron motion varies with the parameters of the problem. Along with 
the amplitude, the importance of cyclotron motion also varies. In the extreme limit where 
the amplitude of the cyclotron motion is very large in comparison with the length scale of 
the problem we study, the whole picture of the dynamical behavior of vortices, as this is 
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determined from the equation of motion for the guiding center alters dramatically. Thus, 
it is important to determine the way the parameters of the model affect cyclotron motion. 
Our numerical investigation revealed a systematic relation between the parameter /3, the 
value of the external current Jq, and the amplitude of the oscillating patterns. In short, 
the cyclotron motion is amplified when (3 decreases as well as when Jq increases. The 
dependence on Jq and [3 is demonstrated in fig.4 and fig. 5 respectively. Note that in the 
runs presented in fig.4 n = 2. and (3 = 1. while in those in fig. 5 k = 1.5 and Jq = 0.03. 
Here we stress once more the analogy between the motion of the vortex in this model and 
the Hall motion of charged particles. The increase of Jq here is equivalent to the increase 
of the external electric field in the charged particle system and the effect on the cyclotron 
motion is the same in the two cases. Also the increase of j3 is physically equivalent to 
the increase of the dielectric constant. The analogue in the particle system is to decrease 
the coupling of the particle with the electric field, without affecting its coupling with the 
magnetic field. Again the consequence is similar in the two systems. 

Finally we consider the case where an initially applied external field is abruptly turned 
off. The response of the vortex to such a "blackout" is of obvious interest. Hall equation 
( |2.2| ) implies that in the absence of external forces the guiding center of the vortex is 
conserved i.e. the vortex is pinned. In other words, while the vortex moves at a constant 
speed just before the current is turned off, it abruptly freezes at the position where it is 
found at the time we switch off the current. Still there is one question to be answered, 
namely, how the location of the MED -which in principle does not coincide with the 
guiding center- will evolve. After all the guiding center is an abstract notion which 
represents the "mean position" of the vortex, while its real position in space is associated 
with the distribution of the energy density. One could possibly assume that after the pause 
of the current the vortex will reorganize itself and will finally relax to a configuration where 
the MED coincides with the pinned guiding center. However, the numerical results lead 
to a completely different picture. 

In fig. 6 (a), (b) we display the results of a simulation of a total duration of 125 time 
units where an external current is on only up to time t = tcrit = 40. The trajectories 
of both the guiding center (solid line) and the MED (dashed line) are plotted. Plot (a) 
shows the trajectories from t = up to t = tcrit while plot (b) displays the trajectories 
till t = 125. As it is shown there, after the current is turned off, the guiding center 
indeed remains fixed at the point where it was found then, while the maximum of the 
energy density sets in a circular motion around the location of the guiding center with a 
radius equal to its distance from the guiding center at the time the current was turned off. 
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The picture described here, has a striking analogy with the Hall effect. A guiding center 
can be also introduced in the case of the two-dimensional electron motion in a uniform 
magnetic field B |^ which again can be interpreted as the "mean position" of the electron. 
When an external electric field is applied the electron sets in a cycloid motion along the 
Hall direction while its guiding center follows a rectilinear orbit along the same direction. 
What is more, when the electric field is turned off, the electron, sets in a circular motion 
around its guiding center which rests. 



5 Discussion 

In this paper, a time-dependent Ginzburg-Landau model for a complex scalar field coupled 
to electromagnetism has been studied mainly numerically. Dissipation was successfully 
incorporated in the model and its effect on the motion of the vortex was analytically 
determined. The results of the numerical study were in accord to the Hall analogy advo- 
cated in ref. based on the derived unambiguous conservation laws. Furthermore, the 
results confirmed with impressive accuracy earlier theoretical predictions concerning the 
speed and the direction of the velocity of the vortex. In short, it was shown that under 
the influence of an external electric current a vortex drifts in a direction opposite to the 
current while in the presence of dissipation it deflects in a direction ranging from 90° to 
180° with respect to the current. 

An important feature which is worth emphasizing is the cyclotron motion i.e. the 
oscillating patterns in the trajectory of the vortex remarkably similar to those encountered 
in the cycloid motion of an electron in the standard Hall effect. Such patterns have been 
already observed in the motion of vortex pairs and in the motion of magnetic bubbles 
in ferromagnetic media |]T^. It seems, that cyclotron motion is a generic feature of 



solitons which exhibit Hall behavior. According to reference there is a whole class 
of field theories whose solitons are expected to exhibit Hall behavior and among them 
there are some interesting variations of the present model W^, |TB[- In some sample 



runs in these models we did encounter cyclotron motion which we consider as a strong 
indication that indeed these systems exhibit Hall behavior. Note that in the framework 
of collective coordinate schemes hke those invoked in |ll|], [|l^ it is not possible to detect 
cycloid patterns, because in the adiabatic limit the amplitude of cyclotron motion becomes 
negligible, a fact which is supported by the results displayed in fig. 4. 

At a phenomenological level, we have presented arguments (mainly by reformulating 
previous results), which are quite encouraging for the relevance of the model to the physics 
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of the superconductor. In particular, we have shown that Hall equation ( p.2|) leads natu- 
rally to the introduction of Magnus force ( |2.5D in the equation of motion for the vortex. 
Furthermore, we argued that the Magnus force has electromagnetic origin due to the in- 
teraction of the moving magnetic flux with the positive ions of the background lattice. 
Finally, we have derived a variation of the Nozieres-Vinen equation for the motion of a 
vortex ( |3.8| ) using plain field theoretical analysis. 

We have also exhibited analytical considerations and numerical results, which suggest 
that in the presence of an external electric current the vortex drifts against the current 
implying a possible link with the opposite sign Hall effect Yet, as we mentioned in 
section IV, the way we incorporate the electric current in the model plays crucial role to 
the derivation of that result. Clearly the next step is to test the predictions of the model 
at hand against more realistic experimental situations. One should find a more natural 
way to introduce the electric current in the specimen. Also one should abandon the 2-dim 
reduction and study the issues presented here in thin films with finite thickness. In such 
a study [^|, preliminary results imply that Hall behavior is also exhibited in the motion 
of magnetic flux tubes which are probed by a surface external current in a 3-dim film 
i.e. a current which is non-zero only at the upper and the lower layers of a 3-dim grid. 
At the experimental front, one might try to mobilize vortices not by applying an electric 
potential on the specimen, but by introducing an electric current in a parallel plane just 
above the specimen. 
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FIGURE CAPTIONS 

Figure 1. A plot of the forces acting upon a vortex which moves with constant velocity 
Yl in the presence of a uniform external current Jq and dissipation. 

Figure 2. The trajectory of the guiding center of the vortex under the influence of an 
external electric current and dissipation. The several lines correspond to different values 
of the friction parameter C^. 

Figure 3. The drift velocity V = {Vi,V2) of the guiding center of the vortex for 
Cd = 4. 

Figure 4. The motion of the vortex as determined by its guiding center (dashed 
line) and the location of the maximum of the energy (solid line) for various values of the 
external current. 

Figure 5. The motion of the vortex as determined by its guiding center (dashed 
line) and the location of the maximum of the energy (solid line) for various values of the 
parameter 

Figure 6. The pinning of the vortex, (a) the trajectory of the vortex while the current 
is on, and (b) the subsequent evolution after the current is turned off. 
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